Draft version January 29, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



PROBING THE DAWN OF GALAXIES AT Z - 9 - 12: 
NEW CONSTRAINTS FROM HUDF12/XDF AND CANDELS DATA ^ 

P. A. Oesch- R. J. BouwENS'*, G. D. Illingworth^ , I. Labbe^, M. Franx^, P. G. van Dokkum^, 
M. Trent:"', M. Stiavelli'', V. Gonzalez^, D. Magee^ 

Draft version January 29, 2013 



o 

(N 

in 

(N 



o 

u 

o 



ABSTRACT 

We present a comprehensive analysis of z > 8 galaxies based on ultra-deep WFC3/IR data. We 
constrain the evolution of the UV luminosity function (LF) and luminosity densities from z ^ 11 
to z ~ 8 by exploiting aU the WFC3/IR data over the Hubble Uhra-Deep Field from the HUDF09 
and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area 
WFC3/IR imaging over GOODS-South. Galaxies are selected based on the Lyman Break Technique 
in three samples centered around 2 ^ 9, z ~ 10 and 2; 11, with seven z ^ 9 galaxy candidates, and 
one each at z ~ 10 and z ^ 11. We confirm a new z ~ 10 candidate (with z — 9.8 ± 0.6) that was not 
convincingly identified in our first z ~ 10 sample. The deeper data over the HUDF confirms all our 
previous z > 7.5 candidates as genuine high-redshift candidates, and extends our samples to higher 
redshift and fainter limits (iJieo ~ 29.8 mag). We perform one of the first estimates of the z '--^ 9 
UV LF and improve our previous constraints at z ~ 10. Extrapolating the lower redshift UV LF 
evolution should have revealed 17 z ~ 9 and 9 z ~ 10 sources, i.e., a factor ^ 3x and 9x larger than 
observed. The inferred star-formation rate density (SFRD) in galaxies above 0.7 M0yr~^ decreases 
by 0.6 ± 0.2 dex from z ~ 8 to z ~ 9, in good agreement with previous estimates. The low number 
of sources found at z > 8 is consistent with a very rapid build-up of galaxies across z ~ 10 to z ~ 8. 
From a combination of all current measurements, we find a best estimate of a factor 10 x decrease 
in the SFRD from z ~ 8 to z ~ 10, following (1 -|- z)'^^-"^^^-^ . Our measurements thus confirm our 
previous finding of an accelerated evolution beyond z ~ 8, and signify a rapid build-up of galaxies with 
Mjjv < —17.7 mag within only ~ 200 Myr from z 10 to z ^ 8, in the heart of cosmic reionization. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function 



> ; 1. INTRODUCTION 

(N ; The launch of the WFC3/IR camera in 2009 signi- 
, fied a major milestone in our ability to observe galax- 
' ' ies within the cosmic reionization epoch at z > 6. 

Thanks to its ~ 40 times higher efficiency for detecting 
T-H , galaxies in the near-infrared (NIR) compared to previ- 
• ous cameras on the Hubble Space Telescope (HST) we 
' have pushed the observational frontier to within only 
~ 450 Myr from the Big Bang. In its first year of op- 
^ . eration WFC3/IR resulted in t he detection of ^ 13 
■ new galaxies at z > 6 (see e.g. iBouwens et al.|[2011b[ ). 
^ ; Three years of science operations of WFC3/IR and sev- 
$— I ■ eral deep extra-galactic surveys have now resulted in a 
large sample of more than 200 galaxies in the reion- 
ization epoch, prima rily at z ^ 7 and z ~ 8 (e.g. 
Bouwens et alJl2011bl: lOesch et a l."2012b' : iMcLure et al.l 



2012t ISchenker et al.l l2012al: iLorenzoni et al.l 120131: 
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i Bradlev et all 12 012: 'Y an et al.l 120111 : iFinkelstein et al.l 
l20ia iGrazi"t al...2012l ) 

From these samples it has become clear that the build- 
up of galaxies during the first Gyr was a gradual process 
at z < 8. The end of reionization left no noticeable 
imprint on the galaxy population (at least down to the 
current limits of Muv — —18 - corresponding to a star- 
formation rate of SFR ~ lM0yr~^). The build-up of the 
UV luminosity function (LF) progresses smoothly across 
the z ~ 6 reionization boundary, following a constant 
trend all the way from z 8 to z ^ 4. Galaxies typically 
become brighter by ^ 30 — 40% per unit redshift accom- 
panied by a proportional (but somewhat larger) increase 
in the average star-formation rate of galaxies (see e.g. 
ISmit et al.ll20Tl IPapovich et al.ll20Tl . 

Given the large samples of galaxies discovered dX z ^ 
7 — 8, the current observational frontier is at z ~ 9 
and at earlier times. This is a period when significant 
evol ution of the galaxy UV LF is e xpecte d from models 
(e.g.lTrenti et aljboiol: iLacev et aLll2011l: iFinlator et al.l 
1201 1| ). The observational evidence has been suggestive 
of a significant drop on the UV luminosi ty density (LD) 
at z > 8, but has not be en conclusive (jBouwens et al.l 
I2011at lOesch et al.ll2012al ). As a result, the extent to 
which the luminosity function and the star formation 
rate density are evolving at z > 8 has been the subject 
of so me debate (see, e.g., iCoe et aI1l2013t iZheng et al.l 
IMl) . 

At these early epochs, current galaxy samples are still 
very small as HST is approaching its limits. Initially, 
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only one z '-^ 10 galaxy candidate was identified (UDFj- 
39546284), even in extremely deep WFC3/IR imaging 
of the Hubble Ult ra Deep Field (HUDF ) as part of the 
HUDF09 survey ()Bouwens et all 1201 lal ). When com- 
bined with all the existing data over the Chandra deep 
field South (CDFS), this one source suggested that the 
galaxy population is changing quickly, building up very 
rapidly from z ~ 10 to z ~ 8. In galaxies with SFR 
> lM0yr~^ (equivalent to Muv — —18 mag), the in- 
ferred UV luminosity density (LD) was found to increase 
by more than an order o f magnitude in only ~ 200 Myr 
from z ~ 10 to z ~ 8 (jOesch et al.ll2012al ). This is a 
factor ~ 5 times larger than what would have been ex- 
pected from a simple extrapolation of the lower redshift 
trends of the UV LF evolution to z ~ 10. 

Several datasets have allowed us to improve these 
first constraints. The multi-cycle treasury program 
CLAS H (PI: Postman) has provided four sources at 
z^9 (|Zheng et al.ll20ia iCoe et aLllMl iBouwens et"all 
l2012al) . In particular, the detecti ons of three z ^ 9 galax- 
ies around CLASH clusters by IBouwens et al.l ()2012aD 
have provided a valuable estimate of the luminosity den- 
sity in this key redshift range. Estimating volume den- 
sities from highly-magnified sources found behind strong 
lensing clusters is challenging, involving systematic un- 
certainties due to the lensing model. IBouwens et al.l 
(|2012af) used a novel technique of comparing the z 9 
source counts to those at z 8 in the same clusters 
and obtained a good relative luminosity density estimate. 
Sinc e the z ~ 8 density is well-established fro m the field 
(e.g. lOesch et al.ll2012bl: iBradlev et aLll2012D . this gave 
a more robust measure than trying to infer source den- 
sities directly using lensing models. Interestingly the 
three z ~ 9 candidates from'Bouwe ns et al.l (|2012al ) are 
completely consistent with the observed drop in the UV 
luminosity density and an accelerated evolution of the 
galaxy population that w as previously seen at z > 8.5 in 
the HUDF and CDFS bv lOesch et all (|2012aD . 

Two other high redshift sources detected in the CLASH 
dataset (one of which is in co mmon with and proceeded 
the sample of Bou wens et a l. 2012a) have added to the 
availa ble constraints. lCoe~t al. (2013) and iZheng et al.l 
(|2012| ) discovered two highly magnified z ~ 10 (z 9.6 
and z ~ 10.7) galaxies in the analysis of the CLASH clus- 
ter data. The luminosity densities inferred from these 
galaxies are somewhat higher, but are very uncertain. 
The large errors on the luminosity density from these 
two detections encompass a wide range of possible trends 
from z ~ 8 to earlier times. However, as we show later in 
this paper, taken together, the sources from the CLASH 
dataset along with the latest sources and constraints 
from the HUDF/CDFS region, are consistent with our 
earlier estimates of substantial accelerated change in the 
luminosity density from Z'^10toz~8. 

Additional progress in exploring the galaxy popula- 
tion at z > 8 has been made through gamma-ray burst 
(GRB) afterglow observations. The current record holder 
of an independently confirmed redshift measurement was 
achieved at z ~ 8.2 fo r GRB090423 (Tanvir et al. 20Q2; 
ISalvaterra et al.ir2009f) . and GRB reds hifts were photo- 
metr ically measured out to z = 9.4 (jCucchiara et al.l 
1201 If ). These measurements can provide additional 
constraints on the total star-formation rate density in 
the very early universe, since GRB rates are thought 



to be an unbias ed tracer of the total star-formation 
rate density fe.g. iKistler et al.l l2009t iTrenti et al.1 120121: 



iRobertson fc EUii 



1201 21) . 



While the initial results are encouraging, it is clear that 
our understanding of the galaxy population at z > 8 is 
still far from complete. Given the very small number of 
sources in each study, it is perhaps not surprising that the 
UV luminosity density measurements at z > 8 are cur- 
rently all within 1 — 2(7 of each other. A next step forward 
in exploring the z > 9 universe can now be taken thanks 
to the 128 orbit HUDF12 campaign (PI: Elhs, G012498). 
While the critical i?i6o observations to discover z ~ 10 
galaxies only reach deeper by ~0.2 mag compared to the 
previous HUDF09 image, the HUDF 12 survey adds deep 
F140W (Ji?i4o) imaging. This allows for Lyman Break 
Gal axy sample selections at z ^ 9 and z '-^ 1 1 — 12 (see 
also IZheng et al.l 120121 : IBouwens et al.ll2012at iCoe et al.l 
[20T3f L 

In a first analy sis of their proprietary HUDF12 data, 
lEUis et all ()2013D compiled a sample of six extremely 
faint z 8.6 — 9.5 galaxy candidates based on a 
photometric redshift technique. One of these sources 
was already in an ea rlier z ^ 7.2 — 8.8 sample of 
Bouwens et'^ (1201 Ibt ) based on the HUDF09 data set. 



Ellis et al.l ( 20131 ) also re-analyzed our previously de- 
tected z ^ 10 candidate Up F.i-39546284 (jBouwens et all 
I2011al: lOesch et al"] l2012aD . From the three years of 
WFC3/IR if 160 data, it is completely clear now that 
XDFjh-39546284 is a real source as it is significantly 
detected in all three major sets of d ata taken in 2009, 
2010 and 2012 (IBouwens et al.ll2012b[) . Very surprisingly, 
however, UDFj-39546284 appears not to be detected in 
the new F140W image, indicating that this source is ei- 
ther a very extreme emission line galaxy at z ~ 2 or that 
it lies at z ~ 12 with the spectral break of the galaxy at 
^1.6/im (see also lBrammer et al.ll2013l ). 

With all the HUDF12 data publicly available, we can 
now extend our search for z > 9 galaxies to even deeper 
limits and to higher redshifts than previously possible. In 
this paper, we perform a search for z ~ 9 — 11 galaxies 
over the HUDF based on the Lyman Break technique. 
This makes use of the fact that the hydrogen gas in 
the universe is essentially neutral at z > 6, which re- 
sults in near-complete absorption of rest-frame UV pho- 
tons short-ward of the redshiftcd Lya line. Star-forming 
galaxies at z > 6 can therefore be selected as blue contin- 
uum sources which effectively disappear in shorter wave- 
length filters. In Section 13.11 we outline our reasons to 
use a Lyman Break selection instead of photometric red- 
shift selection, as is frequently adopted to id entify very 
high-redshift galaxies in the literature, e.g., in lEllis et al.l 
(2013). 

This paper is organized as follows: we start by de- 
scribing the data used for this study in section [2] and 
define our source selection criteria in section [31 where 
we also present our z > 9 galaxy candidates. These are 
subsequently used to constrain the evolution of the UV 
luminous galaxy population out to z ^ 11 in section HI 
where we present our results. In section [SJ we summarize 
and discuss further possible progress in this field before 
JWST. 

Throughout this paper, we will refer to the HST fil- 
ters F435W, F606W, F775W, F850LP, F098M, F105W, 
F125W, F140W, F160W as B435, ^eoe, «775, ^sso, 1^98, 
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TABLE 1 

The 5ct Depths^ of the Observational Data Used in this Analysis 



Field 


Area [arcmin^] 


B435 




«775 


-^814 


2^850 


^05 


^125 


JHi40 


^^160 


HUDF12/XDF'''<= 


4.7 


29.8 


30.3 


30.4 


29.1 


29.4 


29.7 


29.7 


29.7* 


29.8 


HUDF09-1 


4.7 




29.5 


29.3 




29.3 


29.0 


29.2 




29.0 


HUDF09-2 


4.7 


29.5 


29.9 


29.5 




29.2 


29.0 


29.2 




29.3 


ERS 


41.3 


28.4 


28.7 


28.2 


28.5 


28.0 


27.8<* 


28.2 




28.0 


GOODSS-Dccp'^ 


63.1 


28.4 


28.7 


28.2 


29.0 


28.1 


28.3 


28.5 




28.3 


GOODSS-Wide 


41.9 


28.4 


28.7 


28.2 


28.5 


28.0 


27.5 


27.7 




27.5 



Measured in circular apertures of 0'.'25 radius. 

Improved data relative to lOesch et al.l lj2012al ) for z ~ 10 galaxy search. 

A new version of the op tical data is used here (the XDF data set; lUingworth et al. in prep.) compared 
to lBeckwith et aLl 120061) . which results in an improvement in depth of ~ 0.1 — 0.2 mag. 

The ERS field was imaged with Yqqs rather than with I105 . 
* Only the HUDF12/XDF field includes deep JH140 imaging, which wc require for 2 ~ 11 searches and 
which significantly improves z ~ 9 searches. 




Fig. 1.— The WFC3/IR fields over the GOODS South area used 
in this analysis. The HUDF12/XDF field (dark red) contains the 
deepest optical and NIR data to date, which reach to ~ 30 AB 
mag in several bands. The parallel fields HUDF09-1 and HUDF09- 
2 (also dark red) are only 0.5 — 0.8 mag shallower. The wider 
area data covering the whole GOODS-S field are from the ERS 
(yellow) and the CANDELS programs (orange). All these fields 
include imaging in Y105 (or yo98 over the ERS), as well as J125 and 
^^^160 1 which makes it possible to search for z ~ 10 galaxies. The 
HUDF12/XDF field is additionally covered by very deep J-ffi4o 
imaging, which wc exploit to select z ~ 9 Lyman Break galaxies 
and obtain some of the first limits on the galaxy population at 
z ~ 11. 

^105, -^125, JHuo, HiQo, respectively. We adopt ^Im = 
0.3, = 0.7, H"o = 70 kms-^Mpc- \ i.e. h = 0.7 Mag - 
nitudes are given in the AB system (|Oke fc Gunnl[l983l ). 

2. THE DATA 

The core dataset of this paper is the combination 
of ultra-deep ACS and WFC3/IR imaging over the 
HUDF09/HUDF12/XDF field. We enhance this deep 
dataset by using WFC3/IR and ACS data over both 
HUDF09 parallel fields, as well as all CANDELS and 



ERS data over the GOODS-South field (see Figure 
[T|). These datasets provide valuable constraints on the 
more luminous sources, particularly by providing lim- 
its over a larger area than is covered by the small 
HUDF09/HUDF12 field. Ah these datasets include J125 
and HiQo imaging in addition to deep, multi-band opti- 
cal ACS data, which allows for reliable z ^ 10 galaxy 
selections (see Section [373]) . 

All the WFC3/IR and ACS data are reduced follow- 
ing standard procedures. We subtract a super median 
image to improve the image flatness and we register to 
the GOODS ACS frames. For WFC3/IR data we mask 
pixels affected by persistence using the maps provided 
by STScI. The ACS data are corrected for charge trans- 
fer losses when necessary using the public code provided 
by STScI. All images are drizzled to a final pixel scale of 
0'.'06, and the RMS maps are rescaled to match the actual 
fiux fluctuations in circular apertures of 0'.'35 diameter, 
dropped down on empty sky positions in the images. The 
spatial resolution of the data is ^ 0'.'09 and ~ 0'.'16 for 
the ACS and WFC3/IR data, respectively. 

The individual datasets used for our analysis are de- 
scribed in more detail in the following sections. They 
are furthermore summarized in Table [T] and are shown in 
Figure [TJ 

2.1. HUDF12/XDF Data 

Th e Hubble Uhra-Deep Field (HUDF; iBeckwith et"all 
[2OOI) was imaged with WFC3/IR as part of two 
large HST programs now. The HUDF09 (PI: Ilhng- 
worth: iBouwens et al.ll201l"bl ) provided one pointing (4.7 
arcmin^) of deep imaging in the three filters Yios (24 
orbits), J125 (34 orbits), and i?i6o (53 orbits). These 
data were extended re c ently with th e HUDF12 cam- 
paign (jEllis et al.l 120131 iKoekemoer et al.il2012l ). which 
imaged the HUDF further in Y105 (72 orbits) and i/ieo 
(26 orbits), and additionally added a deep exposure in 
JHi4o (30 orbits). These are the deepest NIR images 
ever taken, resulting in a final 5a depth of ifieo ^ 29.8 
mag (see also Table [ij. 

Since the acquisition of the original optical HUDF ACS 
data, several programs have added deeper ACS coverage 
to this region, mainly as part of parallel imaging. We 
combined all the available ACS data over the HUDF, 
which allows us to improve the backgrounds and also 
to push photometry limits deeper by ~ 0.1 — 0.2 mag. 
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These data, along with the matched WFC3 /IR data from 
aU programs, wih be released publicly as the eXtreme 
Deep Field (XDF) and are discussed in more detail in 
lUingworth et al. (in preparation). 

For longer wavelength constraints we also include the 
uhra-deep Spitzer/IRAC data from the 262 h lUDF pro- 
gram (PI: Labbe; see also lLabbe et al.ll2012|) . which reach 
to ~ 27 mag AB (5cr total) in both [3.6] and [4.5] chan- 
nels. These data are extremely important for eliminating 
lower redshift contaminating sources, particularly inter- 
med iate re dshift dusty and/or evolved galaxies (see sec- 
tion [SXl]). 

2.2. HUDF09 Parallel Fields and GOODS-South 

In addition to the HUDF data, we also include the 
two additional deep parallel fields from the IIUDF09 pro- 
gram, as weh as ah the WFC3/IR data over the GOODS 
South field. The latter were taken as p art of the Early 
Release Science (ERS) program (Windh orst et al.|[2011l) 
and the multi-cy cle treasury c a mpaig i i CANDELS (PI: 
Fabe r /Ferguson; iGrogin et al.l 120111 : iKoekemoer et al.l 
I20T1I ). These data were alread y used for a z ~ 10 s earch 
in our previous analysis from lOesch et al.l (|2012a[ ) . We 
therefore refer the reader to that paper for a more de- 
tailed discussion. However, since our previous analysis 
the acquisition of an additional 4 epochs of CANDELS 
DEEP data was completed, resulting in deeper data by 
about 0.2 mag. These are included now in this paper, 
which will allow us to further tighten our constraints on 
the z ~ 10 LF. 

In the optical, we make use of all ACS data taken over 
the GOODS South field, which includes additional imag- 
ing from supernova follow-up programs. These images 
reach ~ 0.1 — 0.3 mag deeper than the v2.0 reductions 
of GOODS, in particular in the zgso-band. We also re- 
duce and include all the /814 data, which was taken over 
this field. By combining all these datasets we have pro- 
duced what is the deepest optical image to date over 
the GOODS-S field. Such deep optical data is very im- 
portant for excluding lower redshift interlopers in LBG 
samples. 

For constraints from Spitzer/IRAC, we use the public 
data from the GOODS campaign. These exposures are 
23 h deep and reach to ~26 mag (M. Dickinson et al. in 
prep). All these fields are also outlined in Figure [TJ 

3. SAMPLE SELECTION 

Source catalogs are obtained with SExtractor 
(|Bertin fc Arnouta I1996D , which is run in dual image 
mode with a specific detection image, depending on the 
galaxy sample we are interested in. For a ll samples at 
z < 10.5, we use a detection image (jSzalav et al.l 
Il999| ) based on the JHuo and i?i6o bands. For z > 10 
JiJi40-dropout selections we use the i?i6o-band for 
source detection. 

All images are matched to the same PSF when per- 
forming photometry measurements. Colors are based on 
small Kron apertures (Kron factor 1.2), typically 0'.'2 
radius, while magnitudes are derived from large aper- 
tures using the standard Kron factor of 2.5, typically 
0'.'4 radius. An additional correction to total fluxes is 
performed based on the encircled flux measurements of 
stars in the iJieo band to account for flux loss in the PSF 



Lya Redshift 
7 8 9 10 11 12 13 14 15 




Fig. 2.— The HST filter set over the HUDF12/XDF and the 
CANDELS fields together with a representative galaxy spectrum 
ai z = 9.5. Due to the high neutral fraction in the IGM, essentially 
all photons shortward of the redshifted Lya line are absorbed for 
galaxies at z > 6. This effect is used to select such high redshift 
sources based on broad-band photometry. The redshift of Lyo is 
indicated on the top axis. The vertical black dotted lines indicate 
the location of the break at ^ = 8.5, 9.5, and 10.5. 



wings. This correction depends on the size of individual 
galaxies and is typically ^0.2 mag. 

3.1. Advantages of Lyman- Break over Photometric 
Redshift Selections 

In this paper, we adopt a Lyman Break galaxy selection 
to identify galaxies at z > 8.5. The major advantages of 
this approach over a p hotometric r edshif t selected sam- 
ple as is used, e.g., in lEUis et al.l ()2013D are simphcity 
and robustness. The Lyman Break technique provides a 
straightforward and robust selec tion, which is easily re- 
producible by other teams (e.g. ISchenker et a l. 2012b). 
Furthermore, the simplicity of the Lyman Break color- 
color criteria also allows for a straightforward estimate 
of the selection volumes based on simulations. 

In contrast, the photometric redshift likelihood func- 
tions are heavily dependent on the assumed template set 
and even on the specific photometric redshift code that 
is used. Additionally, the photometric redshift likelihood 
functions depend on largely unknown priors which are 
needed to account for the number density of intermedi- 
ate redshift passive or dusty galaxies. A particular prob- 
lem is that most high-redshift photometric analyses give 
equal weight to all templates at all redshift (i.e. they 
adopt a flat prior). This includes faint galaxies with ex- 
treme dust extinction at intermediate redshift or passive 
sources at z > 5, which are unlikely to be very abundant 
in reality. 

A further uncertainty, in particular for high redshift 
sources, is how undetected fluxes are treated in the fitting 
process. This can have significant influence on the lower 
redshift likelihood estimates. 

Given all these advantages, we will therefore select 
high-redshift galaxies using the Lyman Break technique 
and we will determine their photometric redshifts a 
posteriori using standard template fitting on this pre- 
selected sample of LBGs. 

3.2. z ^ 9 Lyman Break Selection 

The addition of deep F140W imaging data over the 
HUDF gives us the ability to select new samples of z 9 
galaxies over that field. As can be seen in Figure [2l 
the absorption due to the inter-galactic neutral hydrogen 
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Fig. 3. — The different color-color selections of our LBG samples. For the 2 ~ 9 sample, we average the Y105 and J125 flux measurements 
to provide for a clean dropout selection which separates low redshift galaxies (dashed yellow to red lines) f rom star-formi ng sources at high 
redshift (blue). The different SF tracks assume a dust extinction of E(B — V) = 0, 0.15, 0.30 mag using a lCalzetti et al.l l(2000!) reddening. 
In all panels, the color selections are indicated as light gray regions. The small black dots (in upper panels) show the full galaxy sample, 
while the large filled black circles represent the high-2 candidate sample. Small magenta circles in the upper panels represent galaxies with 
Xopt < ^-S. As can be seen from their color distribution, the most likely sources to scatter into the z ~ 9 selection are blue sources at just 
somewhat lower redshift (2 ~ 7 — 8) rather than intermediate redshift passive galaxies (yellow line at J_f/i4o — -ffieo > 0.2). The upper 
right panel illustrates the Xopt criterion for the 2 ~ 9 selection, which guards our sample against 2 < 7 sources. Sources with Xopt < 0-0^ 
are limited at that value. As can be seen, the 2 ~ 9 candidates lie in a quite unique region in the color-x^pt plot, with only a few sources 
just outside our selection window. For the higher redshift samples (lower left and right panel), the primary selection criterion is the red 
color in the WFC3/IR filters. However, we additionally check for strong detections in both IRAC bands ([3.6] and [4.5]). The limits on 
colors arc Icr. 



shifts in between the Y105 and J125 fihers at z > 9. For a 
robust Lyman Break selection, we thus combine the Y105 
and J125 filter fluxes in which galaxies start to disappear 
at z 9. Our adopted selection criteria are: 

(^105 + Jl25)/2 - Jffl40 > 0.75 (1) 

{Yio5 + Ji25)/2 - JHi4o > 0.75 + 1.3 x {JH^o - i?i6o) 



S/N{B435 to Z85o) < 2 A xlpt < 2.8 

These criteria (shown in Figure [3]) are chosen to select 
sources at z ~ 8.5 — 9.5. We additionally use a ( J125 — 
H160) < 1.2 criterion to cleanly distinguish our z ^ 9 
and z ~ 10 samples (see next section). 

We only include sources which are significantly de- 
tected in the -ffieo and JHno images with at least 3cr 
in each filter and with 3.5(7 in at least one of the two. 
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Fig. 4. — 3"x3"images of the 2 > 8 galaxy candidates. From left to right, the images show, a stack of all optical bands, Kios, Ji?i4o, 
■.^125, fl^ieOi IRAC [3.6], and [4.5]. The stamps are sorted by dropout sample. The approximate photometric redshift of each source is 
shown in the lower left corner of the optical stacked stamp (see also Table [Sj . 
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Fig. 5. — Spectral energy distribution fits to the fluxes of the 2 ~ 9 yj-dropout candidate galaxies in our sample. The magnitude 
measurements and upper limits (Icr) are shown in red. These also include self-consistent flux measurements in IRAC [3.6] and [4.5]. The 
uncertainties we derive for the IRAC fluxes account for the uncertainties involved in removing contamination from neighboring sources, 
resulting in some variation in the effective depth of the IRAC data. The best-flt SEDs (blue) as well as the best low redshift solution (gray) 
are shown as solid lines. The SED magnitudes are indicated as flUed circles. With the exception of source XDFyj-39446317 all sources have 
best-flt redshifts at 2 ~ 8 — 9.5. The low redshift solutions are evolved galaxies at z ~ 2, for which the Lyman Break is confused with the 
Balmer/4000A break. As is evident, almost all other sources show a non-negligible secondary peak in their redshift likelihood distribution 
around z ~ 2. Due to the best-fltting low redshift solution of XDFyj-39446317, we do not include this source in our analysis of the UV 
LF and SFRD . One such source of contamination is completely consistent with our expectation from photometric scatter simulations (see 
Section |3X2)| . 
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As a cross-check we selected sources based on an inverse- 
variance weighted combination of the J125, JH140, and 
7?i6o images at 5cr. Both selections resulted in the same 
final sample of high-redshift sources, i.e. all selected can- 
didates are > 5a detections. 

In addition to the color selections, we require sources to 
be undetected at shorter wavelengths. In particular, we 
use 2(7 non-detection criteria in all optical bands individ- 
ually. Furthermore, we adopt an optical pseudo-x^ con- 
straint. This is defined as xlpt = Z^i SGN(/i)(/i/cr,;)^. 
The summation runs over all the optical filter bands 
S435, Veo6, «775, 4i4, and zsso, and SGN is the sign func- 
tion, i.e. SGN(a;) = -1 if x < and SGN(x) = 1 if x > 0. 
This measure allows us to make full use of all informa- 
tion in the optical data. We only consider galaxies with 
Xopt < 2.8. This cut reduces the contamination rate by 
a factor ~ 3x, while it only reduc es the selection vo l- 
ume of real sources by 20% (see also lOesch et"alll2012b[ ). 
This is a powerful tool for providing source lists with low 
contamination rates (see also Section 13.5.21) . 

These selection criteria result in seven z ^ 9 galaxy 
candidates in the IIUDF12/XDF dataset. These sources 
are listed in Table [5] and their images are shown in 
Figure IH In Figure El we additionally show the spec- 
tral energy distribution (SED) fits and redshift likeli- 
hood functions for these sources. For comparison, these 
are derived fr om two photometric reds hift cod es, ZE- 
BRA (Fcldnian^t al.|[2006 : Oc sch ct al., i2010bD as well 
as EAZY (jBrammer et all |2008( ). As is evident, the 
vast majority of sources does show a prominent peak at 
2^8 — 9 together with a secondary, lower likelihood peak 
at 2; ~ 2. The best- fit photometric redshifts of these can- 
didates range between z = 8.1 — 9.0, with the exception 
of one source (XDFyj-39446317), which has a ZEBRA 
photometric redshift of only Zphot — 2.2. However, using 
the EAZY code and template set the best-fit redshift is 
found at 2; = 8.6. The photometric redshift likelihood 
function for this source is very wide using both codes. 

From photometric scatter simulations (see section 
13.5.21) . we expect to find 0.9 — 1.1 low redshift contami- 
nants in our z ~ 9 sample due to photometric uncertain- 
ties. Therefore, finding a LBG candidate with such a 
low redshift is not necessarily unexpected. We will thus 
list it as possible candidate in Table [S] However, we will 
exclude XDFyj-39446317 from our determination of the 
UV LF at z - 9. 

3.3. The z ~ 10 Lyman Break Selection 

Galaxies at redshifts approachi ng z ^ 10 start to disap- 
pea r in the J125 f ilter. F ollowing iBouwens et al.l (|2011aD 
and lOesch et al.l ()2012aD . we select z 10 galaxies based 
on very red J125 — ^^leo colors and we use Spitzer/IRAC 
photometry to guard this selection against intermediate 
redshift extremely dusty and evolved galaxies in a sec- 
ond step. This selection process also used JHhq data 
when available (i.e. over the HUDF12/XDF field), and 
was used for all the datasets shown in Figure [T] 

The HST selection criteria are: 

( J125 - i?i6o) > 1.2 A {JHuo - i?i6o) < 1.0 (2) 

S/N{Bi^^ to Y105) < 2 A xlpt < 2.8 

in addition to at least 3(T detections in both -ffieo and 
JHi4o and > 3.5a in one of these. All sources in our 



final list also satisfy a > 5a detection criterion in the 
combined J125 + J-ffi4o + Hieo image. 

The J-ffi4o — HiQo color criterion was introduced to 
distinguish z ^ 10 from z 11 galaxies over the 
HUDF12/XDF field. The other fields, which do not have 
deep J-ffi4o imaging, do not include this criterion. We 
account for this difference in our analysis of the selection 
functions (Section 13. 6p . 

When applying these selection criteria to the 
WFC3/IR+ACS data over GOODS-S, we previously 
identified 16 galaxies which satisfied these criteria. How- 
ever, these are all extremely bright in the Spitzer/IRAC 
bands and are even detectable in the shallow [5.8] and 
[8.0] channel data over GOODS-S (having Hieo - [5.8] = 
2.4 — 4.0 mag). These sources were therefore excluded 
from our z ~ 10 analysis, as their -ffieo to IRAC colors 
were too red for a genuine z ^ 10 galaxy. These are 
most likely z ^ 2 — 3 galaxies with significa nt extinction 
and po ssibly evolved stellar populations (see lOesch et al.l 
l2012aD . 

Even taking advantage of the deeper WFC3/IR data 
that became avai lable over the GAN DELS-South field 
subsequent to the lOesch et al.l (|2012a[ ) analysis, no new 
credible z ~ 10 source could be found. However, 
our selection revealed three potential sources in the 
HUDF12/XDF data. Unfortunately, two of these are 
very close to a bright, clumpy foreground galaxy. Their 
photometry is therefore very uncertain, and it is unlikely 
that they are real high-redshift sources. We nevertheless 
list these as potential sources in Table [5l However, we 
will not use them in the subsequent analysis. 

This leaves us with only one likely z ~ 10 galaxy can- 
didate in all the fields we have analyzed here. This is 
XDFj-38126243, which we had previously identified in 
the first-ep och data of the HUDF 09 as a potential z ^ 10 
candidate (jBouwens et al.ll2011al ). However, it was not 
detected at a significant enough level in the subsequent 
second epoch -ffieo data to indicate at high confidence 
that it was real. As a result it was not included in our 
final sample of z 10 sources from the HUDF09 data. 

The source XDFj-38126243 is now clearly detected 
both in the new -ffieo and in the JHhq data from the 
HUDF12 survey, which clearly confirms its reality. This 
is demonstrated in Figure [SI As can also be seen from 
that figure, the source is extremely compact, consistent 
with being a point source. We can therefore not exclude 
that this source is powered by an AGN, which could also 
explain the possible variability over a timescale of 1 year 
(see lower panel of Fig. 15]). However, the low flux mea- 
surement in the second-year HUDF09 data is still con- 
sistent with expectations from Gaussian noise. Taken 
together, the fiux measurements of all three epochs are 
consistent with the source showing no time variability 
(X' = 2.6). 

The source is not detected in our ultra-deep IRAC 
data, and its colors place it at a photometric redshift 
of 9.8iEJ;^ (see also Tableland Figure [71). 

For completeness, we also list three additional poten- 
tial candidates in the appendix in Table [6l These sources 
lie very close to bright foreground galaxies. Formally, 
they show colors consistent with being at very high red- 
shift. However, they are significantly blended with their 
neighbors, such that the fluxes of these sources can not 
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Fig. 6.— Top - Stamps (3" by 3") of tfie z ~ 10 candidate 
XDFj-38126243 — wfiicfi is likely the highest-redshift source over 
the HUDF12/XDF. The top two stamps are the JH140 and Hieo 
observations. The -ffieo stack is also split in observations at three 
different epochs, shown in the bottom three stamps. 'Epoch 1' 
corresponds to the HUDF09 year 1 data (28 orbits), 'Epoch 2' 
are the HUDF09 year 2 data (25 orbits), and 'Epoch 3' are the 
remaining 31 orbits from the HUDF12 and CANDELS observa- 
tions. The S/N in each band is listed in the lower left. This source 
was initially sel ected as high-redshift c andidate after the first year 
HUDF09 data IBouwens et al.ll2011al) . However, as can be seen, 
the galaxy was only very weakly detected (1.4(t) in the HUDF09 
year 2 data. Nevertheless, the source is clearly visible at > 3.5(t 
in all other epochs, as well as in the JH140 data (only taken from 
the HUDF12 program). The source is therefore clearly real. The 
lower signal detection in the 'Epoch 2' data is still consistent with 
the expectation from a Gaussian noise distribution. Bottom - The 
Hieo magnitude measurement for the three different epochs. The 
number of orbits going into each image is indicated close to each 
datapoint. The fluxes are measured in a circular aperture of 0'.'35 
diameter. The magnitude from the total 84-orbit -ffieo image is 
indicated by the gray line, with errorbars represented by the filled 
gray area. The flux measurements are consistent with no variabil- 
ity in this source (x^ = 2.6). However, an AGN contribution to its 
UV fiux can not be excluded. 

be accurately measured with SExtractor, without a more 
sophisticated neighbor subtraction technique. Addition- 
ally, the close proximity to very bright foreground sources 
casts significant doubt on the reality of these sources, and 
we therefore do not include these in our analysis. 

3.4. z ~ 10.7 Lyman Break Selection 

The addition of JH140 imaging from the HUDF12 data 
also allows for the selection of z > 10.5 galaxies based on 
a red Ji?i4o — ifieo color, as the IGM absorption shifts 
through the center of the Ji?i4o filter (see Figure [3]). We 
therefore search for galaxies in the HUDF12/XDF data 
satisfying the following: 
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Fig. 7. — Spectral energy distribution (SED) fit to the fluxes 
of the galaxy XDFj-38126243. The best-fit redshift is found at 
z = 9.8, which is represented by the blue line. The gray SED cor- 
responds to the best low-redshift solution at 2 = 2.5. Open blue 
squares and gray circles show the expected magnitudes of these 
SEDs. Icr upper limits to undetected fluxes are shown as red ar- 
rows. The inset shows the redshift likelihood function as estimated 
with ZEBRA (blue) and EAZY (orange). Both photometric red- 
shift codes consistently find a prominent peak at z = 9.8 and a 
lower-significance peak at z ~ 2.5. 



{JH 



140 



Hieo) > 1.0 



(3) 



S/N{Bi35 to J125) < 2 A xlpt < 2.8 



In order to ensure the reality of sources in this single- 
band detection sample we require > 5a detections in 

^^160- 

Only one source satisfies these criteria: XDFj- 
39546284. This is our previo us highest redshift cari didate 
from the HUDF09 data (see IBouwens et al.l[20TTal ) . Sur- 
prisingly, with JHnQ — HiQo > 2.3 it has an extremely 
red color in these largely overlapping filters. It is by far 
the reddest source in the IIUDF12/XDF data, and its 
rather dramatic color raises some questions. 

At face value, the extreme color, together with the non- 
detection in the optical data and in our deep IRAC imag- 
ing, results in a best-fit redshift of z = 11.8 ± 0.3 for this 
source. However, such a high redshift would imply that 
the source is ~ 20 x brighter than expe cted for z ~ 12 
sources at the same number density (see IBouwens et aD 
2012b). While a strong Lya emission line could reduce 
the remarkably high continuum brightness, the lack of 
Lya seen in galaxies within the reionization epoch at 
z > 6, indicates that the high fraction of neutral hydro- 
gen in the universe at early times absorbs the majority of 
Lyg photons of these galaxies (e.g. Schcnkcr et al. 2012^ 
Pentericci et al.l[2011l : lCaruana et al.ll2012tlBunker et all 
20131 ). Seeing strong Lya emission from a source at 
z ~ 12 during the early phase of reionization is par- 
ticularly unexpected (although perhaps not impossible). 

Alternatively, the source could be a z ~ 2 extreme 
line emitter. An emission line at ^ 1.6/im would have 
to produce the majority of its iJieo flux, which would 
require extreme equivalent w idths. A possible example of 
such a source is presented in lBrammer et afl (|2013f ). For 
more extensiv e discussions of these al tern ative options for 
this source see lBrammer et alj (|2013[ ) and lBouwens et ahl 
f('2012bl) . 
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Given the uncertain nature of this source, we wih treat 
its detection as an upper hmit of < 1 source in the follow- 
ing analysis, and we will only derive upper limits on the 
luminosity and star-formation rate densities at z ~ 10.7 
from this single-source sample. 

3.5. Sources of Sample Contamination 

In the following sections we discuss possible contami- 
nation of our z > 8 LBG samples. 

3.5.1. Dusty and Evolved Galaxies 

As already pointed out in lOesch et al.l ()2012aD , galaxies 
with strong Balmer breaks, or with high dust obscuration 
are a potential source of contamination for z > 9 galaxy 
searches. In particular, in fields with limited depth in the 
WFC3/IR and optical data, such extremely red sources 
can remain undetected shortward of iJieo and can thus 
satisfy the HST selection criteria. Fortunately, with the 
availability of Spitzer/IRAC over all the search fields in 
this study, such sources can readily be excluded from the 
samples based on a j?i6n — [3.6] < 2 color criterion. 

As shown in lOesch et all (|2012a[ l. the CANDELS data 
contains 16 intermediate brightness sources which satisfy 
our z ^ IQ Ji25-dropout selections. However, these could 
all be excluded based on the IRAC constraints. They are 
all found at -ffieo ~ 24—26 mag, which suggests that such 
red, lower redshift galaxies show a peaked LP. Therefore, 
it is expected that they would be much less of a problem 
as contaminants in our fainter samples. 

Given this expectation, it is particularly interesting 
that we actually did not find the lower luminosity coun- 
terparts of such sources in any of the three deep fields, 
even though our IRAC data are sensitive enough thanks 
to the lUDF program. Although the survey volume of 
our deep data is limited, this further suggests that such 
red galaxies are indeed very rare at lower luminosities. It 
remains an open and interesting question as to the nature 
and redshift of these red -ffieo 24 — 26 mag galaxies 
(the redshift is expected to be low, i.e. z 1 — 3, but ex- 
actly over what redshift range they are seen is still quite 
uncertain) . 

3.5.2. Photometric Scatter of Low- z Sources 

After excluding contamination from intermediate red- 
shift, red galaxies, the next most important source of 
contamination is photometric scatter. Photometric scat- 
ter can cause faint, low-redshift sources to have colors 
and magnitudes such that they would be selected in our 
sample. We estimate the magnitude of this effect with 
simulations using real galaxies based on our photometric 
catalogs. 

In particular, we select all sources with -ffieo magni- 
tudes in the range 24 to 25 and we rescale their fiuxes and 
apply the appropriate amount of photometric scatter as 
observed for real sources at fainter luminosities. We then 
apply our selection criteria to these simulated catalogs in 
order to estimate the contamination fraction. This is re- 
peated 5000 times, which results in reliably measured 
contamination fractions. 

As expected, contamination due to photometric scatter 
is most significant at the faint end of our sample. The 
above simulations show that we do not expect to see any 
contaminant at > 1 mag from the detection limit. In the 



HUDP12/XDP z - 9 galaxy sample, we find that 0.9 
contaminants are expected per simulation. Given that 
we find 7 sources, this signifies a ~ 12% contamination 
fraction. Note that this would have been a factor 3x 
higher (2.6 contaminants expected) had we not included 
our optical measurement. Again this shows clearly the 
power of having deep shorter wavelength data, and the 
effectiveness with which data over a range of wavelengths 
can be used. 

Por the higher redshift samples, wc estimate 0.2 and 
< 0.1 contaminants in the HUDP12/XDP LBG selection 
at z 10 and z ~ 10.7, respectively, from analogous sim- 
ulations. Overall, our extensive simulations show that 
the contamination due to photometric scatter is thus ex- 
pected to be < 20% for all these samples. 

As an additional test of the contamination rate in our 
samples, we can use the best-fitting low-redshift SEDs 
for our z 9 candidates in order to estimate with what 
probability such types of galaxies would be selected as 
LBGs. In particular, we use the expected magnitudes 
of the z 2 SED fits shown in Pigure \E\ perturb these 
with the appropriate photometric scatter, and apply our 
z ~ 9 LBG selection. We repeat this simulation 10^ 
times for each of our z ~ 9 galaxy candidates, which 
allows us to estimate the probability for our z ~ 9 sample 
to contain a certain number of contaminants. We find 
that at 65% confidence, our sample contains zero or one 
contaminant, while < 2 contaminants are found in 90% 
of the realizations. Pinally, the chance that the majority 
(i.e., > 3) of these z ~ 9 candidates lie at z 2 is 
estimated to be < 1.5%. 

Note that the average number of contaminants per re- 
alization is found to be 1.1, i.e., very similar to our previ- 
ous estimate of 0.9 contaminants based on using the real, 
bright galaxy population. It should be noted that the 
SED-based test makes no assumptions about the relative 
abundance of faint, star-forming z > 8 galaxies and the 
possible intermediate redshift passive sources at the same 
observed magnitude (~ 29 mag AB). Both observation- 
ally and theoretically, the number density of low-mass, 
passive galaxies at z > 2 is still very poorly understood, 
making it very difficult to gauge the contamination rates 
for our z > 8 samples. In particular, these passive z 2 
galaxies would need to have only 1 — 3 x 10* Mq and 
Mb = -15 to -16 mag. 

Nevertheless, these tests show that while we can not 
completely exclude intermediate redshift contamination 
in our samples, the majority of our candidates are clearly 
expected to lie at z > 8. As we noted earlier, to account 
for this contamination we exclude the z ~ 9 candidate 
XDPyj-39446317 from the subsequent analysis, in agree- 
ment with its ambiguous photometric redshift at z = 2.2 
or z = 8.6. 

3.5.3. Additional Sources of Contamination 

Galactic dwarf stars can be a significant concern for 
z 7 — 8 galaxy samples, due to strong absorption fea- 
tures in their atmospheres, which causes their intrinsic 
colors to overlap with the high-redshift galaxy selection 
criteria. However, this is not as much of a concern at 
z > 9. Stellar spectra are significantly bluer in our se- 
lection colors than high-redshift galaxies. This can be 
seen in Pigure |3l where we plot the location of the stel- 
lar sequence including M, L, and T dwarfs. Stars with 
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Fig. 8. — The redshift selection function for different LBG sam- 
ples. The new selections used in this paper nicely extend the lower 
redshift Y-dropout samples. The new yj-dropout selection has a 
mean redshift z = 9.0, while the J-dropout sample is expected to 
lie at a mean z = 10.0. Although the JHno — Hieo > 1 color is 
only satisfied for 2 > 11 galaxies, the JHi4o-dropout sample ex- 
tends to significantly lower redshift, and peaks only at 2 ^ 10.6. 
This is mainly due to photometric scatter and due to the relatively 
slow change in JH140 — Hieo color from 2 ~ 10 — 11 (see Figure 
[3. 

intrinsically red colors are therefore not expected to be 
a significant contaminant in our samples. The only pos- 
sibility for such stars to contaminate our selection is due 
to photometric scatter, which we implicitly accounted for 
in our photometric simulations in the previous section. 

Additionally, we can exclude contamination by super- 
novae. We verified that all galaxies in our sample are 
detected at statistically-consistent S/N levels in the im- 
ages taken over a time baseline of about Syr as part of 
the HUDF09 and HUDF12 campaigns. 

For a more extensive review of possible contam ina- 
tion in z > 8 sample s, see also iCoe et all (|2013D and 
iBouwens et aD ()2011bl ). 

3.6. Redshift Selection Functions 

The expected redshift distributions of our LBG sam- 
ples are estimated based on extensive simulations of ar- 
tificial galaxies inserted in the real data which are then 
re-select ed in the same manner a s the original sources 
(see also lOesch et al.ll2007ll2012al ). In particular, we es- 
timate the completeness C(m) and selection probabilities 
S{z,m) as a function of iJieo magnitude m and redshift 
z. ^^^^^^^^^^ 

Following IBouwens et al.l (120031 ). we use the profiles 
of z - 4 LBGs from the HUDF and GOODS observa- 
tions as templates for these simulations. The images 
of these sources are scaled to account for the difference 
in angular diameter distance as well as a size scaling of 



(1- 



The latter is motivated by observational trends 



of LBG sizes with r e dshift across z ^ 3 — 8 (see e.g. 
Ferguson et"aLl 12004 IBouwens et ahl 120041: lOesch et al.l 



2010at lOno et al.ll2012ir 



The colors of the simulated galaxy population are cho- 
sen to follow a distribution of UV continuum slopes 
with B = -2.4 ± OA (see e.g . 'Bouwens et al." 120091 
[20T0t iFinkelstein eFaI|[20Tot [Stenwav et al. .200 ^, and 
aremodulated by the IGM absorption model of iMadatj 
(|1995D over a range of redshifts z = 8 to z = 13. 10, 000 
galaxies are simulated for each redshift bin in steps of 
dz = 0.2, which allows for a reliable estimate of the com- 
pleteness and selection probability taking into account 



the dispersion between input and output magnitudes. 

The result of these simulations enables us to compute 
the redshift distribution of galaxies after assuming a LF 
0(A/). 



p{z) = dN/dz 



dm——S{m,z)C{m)(/){M[m,z]) (4) 
dz 



We assume a baseline UV LF evolution with a = —1.73 
and M^z) = -20.29 + 0.33 x (z - 6), consistent with 
the tre nds found across z ~ 4 to z '-^ 8 (|Bouwens et al.l 
l2011b[ ). The normalization is not relevant for the rela- 
tive distributions. For the K-correction in the conversion 
from absolute to observed magnitude, we use a 100 Myr 
old, star- forming template of Bruzual & Chariot (2003). 

The results of this calculation are shown in Figure [H 
It is clear that the redshift selection functions are signif- 
icantly wider than the target redshift range based on the 
simple LBG color tracks shown in Figure [3] The reason 
for this is simply photometric scatter. The mean red- 
shift (and the width) of our samples are 9.0 (±0.5), 10.0 
(±0.5), and 10.7 (±0.6) for the YJ-, J- and JiJ-dropout 
samples, respectively. 

3.7. Comparison to Previous z > 8 Samples 

In the following sections, we compare our LBG samples 
with previous selections over these fields in the literature. 

3.7.1. Bouwens et al. (2011) z > 8.5 Y -dropouts 

In IBouwens et al.l (j2011a[ ). we already identified three 
possible sources at 8.5 z < IQ. These sources were 
identified based on their red Ylos — ^125 colors. With the 
advent of the IIUDF12 data, all these source are con- 
firmed as valid high-redshift candidates. However, they 
are all weakly detected in the Y105 filter, which results in 
a somewhat lower estimate on their redshift. Neverthe- 
less, one of these sources (UDFy-38135540) is included 
in our present z ^ 9 sample. The other two (UDFy- 
37796000 and UDFy-33436598) do have YJ - JH colors 
of ^ 0.4, which are too blue to be included in our sample. 
Their photometric redshifts are 7.8 and 7.7, respectively. 
The photometry of both these sources are also listed in 
the appendix in Table [SI 

As deeper data becomes available it is not unusual to 
find that the photometric redshifts undergo small shifts 
to lower values, also due to the larger nu mber of sources 
at lower redshifts (Muii oz fc Loebl l2008[ ). The original 
bias to higher redshifts results from the larger photo- 
metric scatter in shallower data, resulting in an overes- 
timate of the Lyman break amplitude. Similar biases 
also affect the photometric r edshift samples , e.g.. com- 
pare the redsh i ft est imates of iMcLure et al.l (|201Clf ) with 
iMcLure et all ()2012D . This is a well-known and well- 
understood effect and should be expected to affect all 
redshift estimates derived from photometric data, re- 
gardless of the proced ures used. This effect explains 
the question raised by I Ellis et all (j2013D regarding the 
slightly lower redshift for these sources. 

3.7.2. Bouwens et al. (2011) z ~ 10 Candidates 

In our previous analysis of the full HUDF09 data over 
the HUDF, we already identified the source XDFjh- 
39546284 as a probable high-redshift source. Based on 
those data and on a plausible evolution of the UV LF to 
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UDF12-4106-7304 




^l^- 9-— The flux o f the z ~ 9.5 candidate HUDF12-4106- 
7304 of lEUis et al.l 120131 ) appears to be significantly boosted by a 
diffraction spike. Top - A 15" view around a the bright foreground 
neighbor of the source HUDF12-4106-7304. The source clearly lies 
exactly along the direction of the diffraction spike caused by the 
bright neighbor. Bottom - Sta mps (4" by 4 ) o f the z ~ 9.5 
candidate HUDF12-4106-7304 of lEUis et all (left) next to 

an image of a diffraction spike of a nearby star (center), with the 
rescaled flux of the diffraction spike subtracted from the UDF12- 
4106-7304 image (right). The latter was derived by fitting the 
core of the bright galaxy to the right (West) of HUDF12-4106- 
7304. This was done for both F140W and F160W. The images 
are centered at the same pixel offset from the nearby bright object 
causing the diffraction spike. The bright object lies 5 arcsec to 
the right (W) in each case. The location of the putative z ~ 9.5 
candidate is marked as red circle. It clearly lies extremely close to 
where the peak flux of the diffraction spike is seen. A hint of a 
linear spike is seen in the JHhq image, running across the source 
location. While a source is still seen in the subtracted image, its 
estimated flux is reduced by a factor ~ 2, making this a 2.8(t 
detection only. This very weak detection evidence, together with 
the near-blending of the source with another foreground galaxy, 
strongly suggests that this is not a real source. Such faint higher- 
order diffraction spikes are a well-known pitfall when pushing the 
data to its limits, in particular in ultra-deep imaging data, which 
is mostly taken at the same rotation angle. 

higher redshift, we expected this source to he at z ~ 10.4. 
Surprisingly however, XDFjh-39546284 was not detected 
in the new JHuq data and so its redshift cannot be 
z ~ 10.4 (see also lEllis et"alll2013[ ). Its nature IS now 
unclear. The best- fit z = 11.8 solution is quite unlikely, 
given what we now know about the evolution of the LF at 
redshifts z - 4 - 9. XDFjh-39546284 is ~ 20x brighter 
than expected for a z ^ 12 galaxy at its n umber den- 
sity (see Figure 4 of 'Bouw ens et al1l2012b[ ). Dramatic 
changes to higher luminosity densities at z > 11 are un- 



likely and so this object presents us wit h an interesting 
conund rum. This is discusse d in de tail in lBouwens et al.l 
(|2012bD and'B rammer et al.l ()2013D . 

After the first half of the HUDF09 data was taken over 
the HUDF in the first year of observations, we had iden- 
tified three potential z ^ 10 sour ces (see the Supplemen - 
tary Information/ Appendix A of iBouwens "eFallboilaD . 
These sources were selected as Ji25-dropouts, very sim- 
ilar to the candidates selected in the present analysis. 
However, the three candidates were not detected at suf- 
ficient significance in the second year WFC3/IR i?i6o 
data, which raised the possibility that they were spuri- 
ous detections. 

With the advent of additional iJieo data from the 
HUDF 12 survey, we can now confirm that all these three 
sources are in fact real. They are all significantly de- 
tected in the full i/igo and JHhq data. However, only 
one of these sources is now in our z > 8 galaxy sam- 
ple. Two sources show photometric redshifts of z ^ 8, 
given their very faint detections in the Fios data of 0.5a 
and 2.3(7, respectively. However, we remark that one 
of these two sources may still be at z > 8.5 given the 
tentative nature of its Y105 band detection (i.e. 0.5ct). 
The last source (XDFj-38126243) remains in our new 
z ~ 10 Ji25-dropout sample. For this source, we find 
a photometric redshift of z = 9.8lQ g. With the possible 
exception of the enigmatic z ~ 11.8 redshift candidate 
(XDFjh-39546284), this is therefore the highest redshift 
galaxy candidate in the HUDF12/XDF field. 

In Figure [51 we show the JHhq and iJieo stamps of 
the source XDFj-38126243, including splits of the data 
by epoch. It is clear that the source is real, as it is now 
detected at 5.8cr in -ffieo and at 3.4cr in JHuq. As can 
be seen, the second year HUDF09 data (Epoch 2), only 
contains a weak, though statistically-consistent, signal 
of this source. Given that there were just two epochs 
available at that time and the overall S/N of the source 
was be low our thresho l d, we d id n ot include thi s source 
in the IBouwens et"all (|2011a[ ) and lOesch et alJ (|2012aD 
analyses. 

Also shown in Figure[6]is the best-fit template and pho- 
tometric redshift distribution for this so urce. With both 
photometric redshift cod es ZEBRA (iFeldmann et al. 
20061 lOesch et all l2010bD and EAZY ()Brammer et al. 
200.q) , we find a consistent best- fit photometric redshift 
at z = 9.8 — 9.9 with uncertainties of Az ~ 0.6 (see Ta- 
ble [5]). As expected for such a faint source, the redshift 
likelihood function shows a lower redshift peak around 
z ~ 2.5 (gray SED). The integrated low-redshift (z < 5) 
likelihood is 18%. This is consistent with our estimate 
of lower redshift contamination due to photometric scat- 
ter in our Ji25-dropout sample. Taken together the data 
are consistent with this being a viable and likely z '--^ 10 
candidate. 

Note that lEllis et all (|2013[ ) do not include this source 
in their analysis, but they specifically discuss it. They 
state the source is not significantly detected (< 5a) in 
their summed J125 -I- J-ffi4o + ^^leo image, for which 
we see two main reasons. (1) With a photometric red- 
shift of z ^ 10, this galaxy is largely redshifted out of 
the J125 band, greatly reducing (by ~ 30%) its detec- 
tion significance in a J125 + Ji?i4o + ^^leo stack. Use 
of a JHi4o + Hiao stack is better for such cases and 
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is what we do here. (2) The source is very compact, 
and therefore it is detected at higher significance in the 
small apertures we use here f or S/N mea surements (0'.'35 
diameter) compared to the lEllis et all (2013) analysis 
(0.47 — 0.50"diameter). Moreover, we stress that this 
source is significantly detected in several independent 
sub-sets of the data, and is therefore certainly real (see 
Figure [6]). 

3.7.3. Elhs et al. (2013) and HUDF12 Team Papers 

The HUDF12 team has recently published a sample 
of seven z > 8. 5 galaxy ca,ndida tes identified in the 
HUDF12 data in lEllis et all (|2013[ ). These sources were 
base d on a photometric redshift selection technique (see 
also iMcLure et al.ll2012l ). with a > 5fT detection in the 
•/i25+</-ffi4o+^^i60 summed image. However, as shown 
in Figure [2] galaxies start to disappear in J125 at z > 9, 
which is why we could expe ct to find additio nal sources 
in our sample compared to lEllis et all ()2013f ) . Further- 
more , we use smaller apertures for S/N measurements 
than lEllis et all ()2013l) , which in most cases are more 
optimal for such very small z ^ 9 — 10 sources. This 
results in a small additional gain of ~ 20% in S/N. 

In general, o ur sample is in ve ry good agreement with 
the selection of lEllis et al.l (|2013D . With the exception of 
two, we include all their sources in our z > 8 samples. 
The discrepant ones are UDF12-3895-7114 and UDF12- 
4106-7304, which we discuss below. Their photometry is 
additionally listed in the appendix in Table |6l 

UDF12-3895-7114: This source certainly show colors 
very similar to a z > 8 candidate. However, we measure 
(YJ) — JHi4o — 0.5 ± 0.5, which is bluer than our selec- 
tion color for the z ^ 9 sa mple. Hence it is not included 
m our z ~ 9 sample. While lEllis et al.l (|2013l ) find a best- 
fit photometric redshift of 8.6 ± 0.7, the source is no t 
present in the 'robust' sample of iMcLure et al.l (|2012l ). 
and we find a photometric redshift distribution function 
which is very broad, with a best-fit at z 0.5 (using 
both ZEB RA or EAZY ). This different resuh compared 
to the Elli s et all ()2013| ) redshift estimate may be caused 
by small uncertainties in the photometry measurements 
(given that we also use different apertures). Addition- 
ally, we perform IRAC flux measurements on a source 
by source basis. These include an additional uncertainty 
due to the subtraction o f neighboring sources (see e.g. 
Labbe et a l. 2010', l 2012t ). This is not the case in the 
Ellis et cdJ (2013) analysis, who note that they use con- 
stant upper limits on the IRAC fluxes. 

As we have discussed, photometric redshifts are very 
uncertain for sources this faint and so there is a chance 
that this source is still at z > 8. Nevertheless, our analy- 
sis raises significant doubt about its high-redshift nature. 

UDFI2-4IO6-73O4: The WFC3/IR PSF shows signif- 
icant diffraction spikes, which are caused by the mount of 
the secondary mirror. While typically only seen around 
bright stars, these diffraction spikes are so strong that 
they can also emanate from compact foreground galaxies, 
particularly in the red der WFC 3/IR filters. The source 
UDF12-4106-7304 of Ellis et al.l (2013) is located at the 
edge of such a diffraction spike for both the i?i6o and 
the JHi4o fihers (the only fibers wherein UDF12-4106- 
7304 is significantly detected). This is shown in Figure 
[9l The photometry of this source is clearly significantly 
enhanced by the diffraction spike. The detection signifi- 
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Fig. 10. — Expected number of z ~ 10 candidates per bin of 0.25 
mag in the different fields considered in our analysis, assuming 
that the UV LF evolves steadily from 2 ~ 8 to z ~ 10, consis- 
tent with the well-established trends from 2 ~ 8 to 2 ~ 4. With 
this assumption the HUDF12/XDF alone should have contained 
4.6 candidate galaxies at 2 ~ 10. In our whole survey area, we 
would have expected to see ~ 9 sources now. Given that only one 
candidate galaxy could be identified (shown by the arrow), this 
provides strong, direct evidence that the UV luminosity function 
and LD are evolving rapidly from 2 ~ 10 to 2 ~ 8. 

cance of UDF12-4106-7304 is critically reduced once the 
diffraction spike signal is removed. The profile of the 
bright foreground galaxy is non-trivial to model, but for- 
tunately only its core is relevant for caus ing the diff rac- 
tion spikes. We therefore use galfit (Peng et al.ll200 2') to 
model the center of this source and subtract the diffrac- 
tion spikes that were scaled to match the core flux. Doing 
so results in the flux of UDF12-4106-7304 being reduced 
by a factor 2x, both in Ji?i4o and Hiqq, which makes 
it only a 2.8ct total NIR detection. This is too low to be 
included in a robust sample. 

Additional uncertainty about the reality of this source 
arises due to its different morphology in the Hiqq and 
JHi4Q images. The 'source' also lies close to another faint 
foreground galaxy. It is therefore not clear whether the 
diffraction spike and the neighboring galaxy conspired to 
lead to the detection of this potential candidate. In any 
case, for these reasons, and for the low detection flux, 
the reality of UDF12-4106-7304 remains in question and 
we do not include this source in our analysis. 

4. RESULTS 
4.1. The Abundance 0/ z > 8 Galaxies 

The sample of nine z > 8 galaxy candidates we com- 
piled in the previous sections allows us to make some of 
the first estimates of the z ^ 9 — 11 UV LFs. Although 
limited in area, the HUDF12/XDF data alone provide 
very useful constraints already at z ~ 9 and limits at 
z 11. Additionally, due to the deeper data over the 
HUDF and the CAND ELS GOODS - S field compared to 
our previous analysis in lOesch et al.l ()2012a|) . we are able 
to improve our constraints on the z ~ 10 LF. 

These new constraints at z ~ 9—11 will allow us to test 
whether the galaxy population underw ent accelerated 
evoluti on at z > 8 as previously found in IBouwens et al.l 
(|2011aD and lOesch etall (|2012aD . or whether the UV LF 
trends from lower redshift continue un changed to z > 8 
(the preferred i nterpretation of, e.g., lEllis et all 120131 : 
ICoe et aLllMl IZheng et al.ll20ia) . 

In order to test for such accelerated evolution, we start 
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by estimating the number of galaxies we would have ex- 
pected to see in our z 9 — 11 LBG samples, if the 
lower redshift trends were to hold unchanged at z > 8. 
By comparing the observed number of sources with those 
expected from the extrapolations we derive a direct esti- 
mate of any changes in the evolution of galaxies at z > 8. 

To derive the expected numbers we use our estimates 
of the selection function and completeness measurements 
described in section [3761 This allows us to compute the 
number of sources expected as a function of observed 
magnitude. For an assumed LF (/)(M), this is given by: 



dm 



Am 



dz^S(m,z)C(m)(t>(M\m,z]) (5) 
dz 



For the UV LF evo lu tion, we adopt the rela- 
tions of iBouwens erall ()2011b[ ): i/)* = 1.14 x 
10^'^Mpc~'^mag~^ =const, a — —1.73 =const and 
M * (z) = -20.29 + 0.33 x (z - 6). Note that we as- 
sume constant values for the faint-end slope a and the 
normalization These relations are used as a baseline, 
when extrapolated to higher redshifts, to test whether 
the observed galaxy population at z > 8 is consistent 
with the trends at later times, i.e., at lower redshifts. 

With these assumptions, we find that we would ex- 
pect to detect a total of 17 ± 4 galaxies in our z ~ 9 
yj-dropout sample over the HUDF12/XDF field alone. 
Yet, after correcting for one potential contaminant (see 
Section [3.5. 2|) . we only detect six sources. This is 2.8 x 
fewer than expected from the trends from the lower red- 
shift LFs. 

Similarly, with the same assumptions about the evolu- 
tion of the UV LF from Z'^4toz'^10, we would expect 
to see a total of 9±3 sources in our z ~ 10 J-dropout sam- 
ple. We only find one such source, which suggests that 
beyond z 9, the decrement compared to the baseline 
evolution is even larger than we found previously from 
the HUDF09 and the 6-epoch CANDELS data. We now 
expect to see three more sources compared t o the six that 
we exp ected to see in the earlier analysis of lOesch et al.l 
(|2012aD . Yet, no additional z ~ 10 sources are found. 

The expected magnitude distribution for the z ^ 10 
sample is shown in Figure [TUl As indicated in the fig- 
ure, our search for z ~ 10 sources in the CANDELS 
and ERS fields of GOODS-S should have resuhed in 
two detections. Only 0.5 sources were expected in these 
fields from our previous analy sis using somewhat shal- 
lower data (|Oesch et al.ll2012aD . 

Depending on the assumptions about halo occupation, 
the expected cosmic variance for a sin gle WFC 3/IR field 
is 40 - 45% (Trentlfe StiavolH 2008| [Robertson 2010). 
In order to estimate the significance of our finding of a 
large offset (i.e., decrement) between the expected num- 
ber of sources and that seen, we have to combine the 
Poissonian and cosmic variance uncertainties. We esti- 
mate the chance of finding < 1 source in our full search 
area using the appropriate expected number counts and 
cosmic variance estimates for the individual search fields. 
T he latter are based on th e cosmic variance calculator 
of iTrenti fc Stiavellil (|2008[ ). Using simple Monte-Carlo 
simulations, we derive that given that we expected to 
find 9 sources, finding < 1 occurs at a probability of only 
0.5%. Therefore, our data are inconsistent with a sim- 
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Fig. 11. — Constraints on the 2^9 (top) and 2 ~ 10 (bottom) 
UV LF from the HUDF12/XDF data as well as from the additional 
fields for z ~ 10. Lower redshift LFs are shown as gray solid lines 
for illustration of the LF evolution trends. T hese are the most re- 
cent determinations from lBouwens et al.l I|2007i . l2012a 'l at 2 = 4 — 7 
and Oesch et al. (20121^) at 2 ~ 8. (Top) - Our step -wise 2 '-^ 9 
LF (dark blue circles) is computed in bins of 1 mag, which contain 
1, 2 and 3 sources, respectively. These measurements are consistent 
with (but consistently below) the expected LF given an extrapola- 
tion from lower redshift (dashed blue line). The best-fit LF based 
on luminosity evolution is shown as a solid blue line. This is de- 
rived from the expected number of sources in bins one magnitude 
wide, and is a factor ~ 1.5 — 4x below the extrapolated LF. Also 
shown a s green squares is the step-wise 2 ~ 9 determination from 
IMcLure et al. ( 2013) who use a photometric redshift sample de- 
rived from HUDF12 data. Their determination is in very good 
agreement with our measurement, although it is unclear why they 
only constrain the LF at the very faint end. The dot-dashed line 
represents the best z ~ 9 LF estimate ofLBouwens et al. (2012^) 
over the magnitude range where it is constrained by their Icnsed 
candidates from the CLASH dataset. Their determination is based 
on scaling the normalization of the 2 ~ 8 LF to account for the 
low number density of z ~ 9 LBG candidates found over the first 
19 clusters. Within the current uncertainties, all three determina- 
tions of the 2 ~ 9 LFs are in good agreement, finding accelerated 
evolution compared to the lower redshift trends. (Bottom) - At 
2 ~ 10, our analysis includes several additional fields, which is why 
we can probe to much lower volume densities than for our 2 ~ 9 
LF. Nevertheless, since we only find one potential 2 ~ 10 galaxy 
candidate in our data, we can mostly only infer upper limits on 
the LF. Again, these are consistently below the extrapolated LF 
(dashed red line), indicating that the galaxy population evolves 
more rapidly at 2 > 8 than at lower redshift. 



pie extrapolation of the lower redshift LF evolution at 
99.5%. 

This new estimat e re inforces the concl usion of 
IBouwens etHI (|2011aD and lOesch eFaLl (|2012al) that the 
evolution in the number density of star forming galaxies 
between z ^ 10 and z ~ 8 is large, and larger than ex- 
pected from the rate of increase at later times, i.e., the 
evolution was accelerated in the ^ 200 Myr from z ~ 10 
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TABLE 2 
Stepwise Determination of the 
z 9 AND 2 ~ 10 UV LF Based on 
THE Present Dataset 



Muv [mag] 


<t>t [10~^Mpc-3mag-l] 




2 ~ 9 


-20.66 


< 0.18 


-19.66 




-18.66 


U.OO_Q 24 


-17.66 


1 6+°-^ 


2 ~ 10 



-20.78 < 0.0077* 

-20.28 < 0.013 

-19.78 < 0.027 

-19.28 < 0.083 

-18.78 < 0.17 

-18.28 < 0.34 

-17.78 0-58lo:5o 

1(T upper limit for a non-detection. 

to Z ~ 8. 

4.2. The UV Luminosity Function at z > 8 

The above calculations of the expected number of 
sources can be used directly to constrain the UV LPs 
at z ~ 9 — 11. Since the number of candidates are small, 
we need to make some assumptions about how to charac- 
terize the evolution. The UV LPs at later times provide 
a valuable guide. The parameter that evolves the most is 
the characteristic luminosity L*; the normalization 
and the faint-end slope (a) are relatively unchanged from 
z ~ 8 to z 4 (although we do have evidence for evo- 
lution toward steeper faint-end slopes). This suggests 
that we should estimate what evolution in the charac- 
teristic luminosity best reproduces the observed number 
of sources, while keeping both the normalization and the 
faint-end slope fixed. 

Doing so results in a best-estimate of the luminosity 
evolution of dM/dz = 0.49 ± 0.09 from z ~ 6 to z ~ 9 
and dM/dz = 0.6 ± 0.2 to z ^ 10. The characteristic 
magnitudes at these redshifts are thus expected to be 
Af,(z - 9) = -18.8 ±0.3 and M,(z - 10) = -17.7±0.7. 
The uncertainties on these measurements are still quite 
large, given the small sample sizes and the small area 
probed, in particular for the z 9 search. 

We perform the same calculation for the z ~ 10.7 
JiJi40-dropout sample. However, given the uncertain 
nature of the single candidate source in that sample 
(XDPjh-39546284), we treat the estimate at z - 10.7 
as an upper limit. We therefore compute the evolution 
in which is needed to produce one source or fewer 
in the sample. This is found to be dM/dz > 0.4, which 
is a less stringent constraint than that for our z ^ 10 
estimate, due to the much smaller area probed by our 
z ^ 10.7 search. The inferred constraint on the charac- 
teristic magnitude is M*(z — 10.7) > —18.4 mag. All 
our estimates of the UV LPs are summarized in Table [3l 

The above estimates for the characteristic luminosity of 
the UV LP can also be compared with the characteristic 
luminosity from the step-wise determination of the LP 
using the observed galaxies and limits. The step-wise 



luminosity function is derived using an approximation of 
the effective selection volume as a function of observed 
magnitude Vce{m) = dz^S[z,m)C{m). The LP is 
then simply (j){Mi)dM = N°^yVafiim,). 

This derivation is only valid as long as the absolute 
magnitude varies slowly with observed magnitude. This 
is not the case for z > 11, since the IGM absorption 
affects the -ffieo band such that the relation between lu- 
minosity and the observed magnitude becomes strongly 
redshift dependent. We therefore restrict our analysis 
of the step-wise LPs to the z < 10.7 samples. In any 
case, the small area probed by the HUDP12/XDP data 
does currently not significantly constrain the UV LP at 
z ~ 11. The step- wise z ^ 9 and z '-^ 10 LPs are tabu- 
lated in Table O 

PigurefTTlshows our constraints on the UV LP at z '-^ 9 
and z ~ 10. The expected UV LPs extrapolated from the 
lower-redshift trends are shown as dashed lines. Clearly, 
at both redshifts, the observed LF lies significantly below 
this extrapolation, as expected from our analysis of the 
observed number of sources in the previous section. 

The best-fit z 9 LP using A/* evolution is a factor 
1.5 — 5 X below the extrapolated LP at Muy > —20 mag. 
At the bright end, the small area probed by the single 
HUDP12/XDP field limits our LP constraints to > 10"'' 
mag~^Mpc^"^, which is clearly too high to be meaningful 
at Mjjv < — 20. It will therefore be very important to 
cover several fields with P140W imaging in the future in 
order to push the z ^ 9 selection volumes to interesting 
limits. 

At z ^ 10, the use of deeper data both on the HUDP 
and on the GOODS-S field allows us to push our previous 
constraints on the UV LP to fainter limits. Since we only 
detect one z ~ 10 galaxy candidate, however, our con- 
straints mainly consist of upper limits. Nevertheless, it 
is evident that at all magnitudes Muy > — 20 mag these 
upper limits are consistently below the extrapolated UV 
LP, up to a factor 4x lower. Additionally, the limits are 
also clearly below the best-fit z ~ 9 LP, showing that the 
UV LP continues to decline at z > 9. 

4.3. The UV Luminosity Density Evolution at z > 8 

The evolution of the UV luminosity density (LD) at 
z > 8 has received considerable attention in recent pa- 
pers, triggered by our initial finding of a significant drop 
in the LD from z ~ 8 to z 10 (i.e., a rapid increase 
in the LD within a short period of time). With the new 
IIUDP12/XDP data, it is now possible to refine this mea- 
surement by adding a z '-^ 9 and a z ^ 11 point, while 
also allowing us to improve upon our previous measure- 
ments at z ^ 10. 

The UV LDs inferred from our z > 8 galaxy sam- 
ples are shown in Pigure 1121 The measurements show 
the LD derived by integrating the best-fit UV LP de- 
termined in the previous section. The integration limit 
is set to Muv — —17.7 mag, which is the current limit 
probed by the IIUDP12/XDP data. Por comparison, we 
also show the lower redshift LD measurem ents from the 
compilation of lBouwens et all (|2007l . l2012bD . These were 
computed in the same manner as the new z > 8 values, 
and were not corrected for dust extinction. A summary 
of our measurements for the LD are listed in Table IH 

As can be seen, our new measurements at z > 8 lie 
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TABLE 3 

Summary and Comparison of z > 8.5 LF Determinations in the Literature 



Reference 


Redshift log 


4>t [Mpc ^mag -*-] 


M*^ [mag] 


a 


This Work 


11 


-2.94 (fixed) 


> 18.4 (Icr) 


-1.73 (fixed) 


This Work 


10 


-2.94 (fixed) 


-17.7 ±0.7 


-1.73 (fixed) 


This Work 


9 


-2.94 (fixed) 


-18.8 ±0.3 


-1.73 (fixed) 


Oesch et al. f2012a~) 


10 


-2.96 (fixed) 


-18.0 ±0.5 


-1.74 (fixed) 


Bouwens et al. (2012a) 


9.2 


-3.96 ±0.48 


-20.04 (fixed) 


-2.06 (fixed) 
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Fig. 12. — The evolution of the UV luminosity density (LD) 
PUV contributed by all galaxies brighter than Afi4oo = —17.7 mag. 
Our new measurements from the HUDF12/XDF data are shown 
as red squares. Measu rements of the LD at z < 8 are derived 
from the UV LFs from IBouwens et al.l | [2007l . I2012d) . No correc- 
tion for dust extinction has been applied. The measurements are 
plotted at the mode of the redshift distributions shown in Fig- 
ure 6. For the highest redshift z ~ 10.7 JH-dropout sample, we 
only show an upper limit given that the single source we find in 
this sample is either at an even higher redshift (where the selec- 
tion volume of our data is very small) or is a low redshift extreme 
line emitter. The dark gray line and shaded area represent an ex- 
trapolation of the redshift evolution trends of the z ~ 4 — 8 UV 
LF. Our LD estimates at z ~ 9 — 10 are clearly lower than this 
extrapolation. However, the observed rapid build-up of galaxies 
at z ~ 10 to z 8 is not unexpected, since it is consistent with 
a whole suite of theoretical models. Some of these are shown as 
colored li nes. They are halo occupation models (blue solid and 
dashed, Trenti et al. 2010; Tacchella et al. 2012), a semi-analytical 
model (orange dashes, Lacey et al. 2011), and two hydrodynamical 
simulations (Fin latbr et'al.ll2011l : lJaacks et al.ll20'l2l 'l. These difl'er- 
ent models uniformly predict a steepening in the LD evolution at 
z > 8. The conclusion to be drawn is that the shape of the trend 
from z ~ 10 to z ~ 7 is mainly due to the rapid build-up of the un- 
derlying dark-matter halo mass function, rather than any physical 
changes in the star-formation properties of galaxies. 



significantly below the z ^ 8 value. The decrement in 
LD from z ~ 8 to z ^ 9 is 0.6 ± 0.2 dex, and it is 
even larger at 1.5 ± 0.7 dex to z ^ 10. Therefore, our 
data confirms our previous finding of more than an order 
of magnitude increase of the UV LD in the short time 
period, only 170 Myr, from z ~ 10 to z ~ 8. 

The gray line and shaded area show the expected LD 
evolution when extrapolating the z ~ 4 — 8 Schechter 
function trends to higher redshift. All our measurements 
at z > 8 lie below the extrapolation. Although the offsets 
individually are not large (they are < 2a) , the consistent 
offset to lower LD supports a hypothesis that significant 
changes are occurring in the LD evolution at z > 8. 

It is interesting to note that this offset to lower LD 



is not unexpected, as it is also seen in several theoreti- 
cal models. In Figure [T^ we compare our observational 
resul ts to two conditio nal l uminosity funct i on m odels 
from lTrenti et all (|2010[ ) and lTacchella et al l (I2012D. the 
predic tion from a semi-analvtical model of iLacev et al.l 
(|2011| ). as well as the results from t wo h ydrodynami- 
cal sim ulations of iFinlator et al.l (j2011[) and ljaacks et alj 
(|20T1 . 

All these models are in relatively good agreement with 
the lower redshift (z < 8) measurements. As can be seen 
from the figure, essentially all models do show a steeper 
evolution at z > 8 than a purely empirical extrapola- 
tion of the UV LF further into the epoch of reionization. 
Since these models are all very different in nature, this 
strongly suggests that the rapid build-up we observe in 
the galaxy population is mainly driven by the build-up 
in the underlying dark matter halo mass function, which 
is also evolving very rapidly at these epochs. 

4.4. The SFR Density Build-up from Z'^11 to z ^ 8 

Our results combined with those of others now provide 
a substantially larger sample at z > 8 for estimate of 
the SF R de nsity than w a s av ailable for IBouwens et al.l 
(|2011aD and lOesch eFalj (p0i2a). The SFR density at 
z > 8 was recently estimated based on fo ur high-redshift 
galaxi e s identified in the CLASH survev (IBouwens et al.l 
l20T2i ICoe et all I20TI IZheng et al.l |201J), and from 
seven galaxies identified in the HUDF12 data ([Ellis et al.l 
120131 ). We present all these results, together with our 
own measurements in Figure 1131 where we plot the SFR 
density as a function of redshift in star-forming galaxies 
with SFR > O.7M0yr^^ (corresponding to a magnitude 
limit of Mjjv = —17.7 mag). 

The SFR densities are derived from the UV LD esti- 
mates after correction for dust extinction. We use the 
most recent determinations of the UV continuum slopes 
/3 as a function of U V luminosity and redshift from 
IBouwens et al.1 (|2012dl ) together with the iMeurer et al.l 
19991) /3-extinction relation. The dust-corrected LDs are 



then conv erted to SFR den sities using the conversion fac- 
tor of M adau et al.l (|199g), assuming a Salpeter initial 
mass function. 

Across z ~ 4 to z ~ 8, the star-formation rate density 
clearly evolves very uniformly. The evolution is well re- 
produced by a power law oc (1 -I- 2;)-3-6±o.3^ which is 
shown as dark gray line in Figure [TS] Interestingly, all 
measurements lie below the extrapolation of this trend 
to higher redshift. Again, individual measurements are 
within < 2a of the trend, but the offset to smaller SFR 
densities in the mean is very clear. 

Note that each of the three groups find a consistent 
decrement of the SFRD from z ~ 8 to z '-^ 9. Specifically, 
from our data, we find a drop by 0.6 ± 0.2 dex. This is 
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TABLE 4 

Summary of Luminosity Density and Star-Formation Rate Density 

Estimates 



Dropout Sample 


Rcdshift 




logio P* 






[erg s"lHz~^Mpc^^] 


[MQyr-lMpc~3] 


YJ 


9.0 


25.00«-i? 




J 


10.0 


24.ll°i 


Q 7+0.7 
' -0.9 


JH 


10.7 


< 24.6 


< -3.3 


B 


3.8 


26.38 ± 0.05 


-1.21 ±0.05 


V 


5.0 


26.08 ± 0.06 


-1.54 ±0.06 


i 


5.9 


26.02 ± 0.08 


-1.72 ±0.08 


z 


6.8 


25.88 ± 0.10 


-1.90 ±0.10 


Y 


8.0 


25.58 ± 0.11 


-2.20 ±0.11 



Note. — Th e lower rcdshift data points ar e bas ed on the UV LPs fr om 
IBouwens eFall (|2007), Bouwcns et al.l l l2011bl ). and lOesch et al.l ll2012bl) 



t Integrated down to 0.05Z/*_3 (Afi4oo = —17.7 mag) 



O 

Oh 







Q 

o 



-1 

-1.5 
-2 

-2.5 
-3 

-3.5 
-4 

-4.5. 




> 0.7 M^/yr 



HUDF09/12+GOODSS 

Ellis+13 
CLASH 



4 



6 7 
Redshift 




8 9 10 



Fig. 13. — The evolution of the star-formation rate density (SFRD) p, contributed by star-forming galaxies brighter than M1400 = —17.7 
mag. This limit corresponds to a star-formation limit > O/7M0yr~^. Measurements at ^ > 8 are shown from the present analysis (dark 
red squares), from the HUDF 12 anal ysis of Ellis ct aLl 1120131 ) (green circles) as well as from CLASH cluster detections (blue triangles; 
luwens et al.l 12012a': 'Coe ct al. '201^ [Zheng et al.ll2012j ). The lower redshift datapoints are derived from UV LFs from IBouwens et al.l 
307, 2012c). The SFRD is derived from the UV LD integrated to these limits and corrected for dust extinction using the most recent 
estimates from iBouwens et al.i (i2012dl) . A clear decrement by 0.6 dex in the SFRD is consistently seen between the measurement at 2 ~ 8 
and that at z ~ 9. At redshifts higher than 2 ~ 9, all datapoints only contain one galaxy, resulting in large uncertainties. (We have 
corrected down the measurement of Ellis et al. to account for a likely diffraction spike source; c.f. open vs. filled circle at z ~ 10.) Given 
the large uncertainties, the individual z > 8 measurements are all consistent with each other. At z = 4 — 8 the SFRD increases gradually, 
following p, oc (1 + z)~^'^^'-^-^ (dark gray line). The extrapolation of this trend to higher redshift is shown by the dashed line and gray 
region (Ic). All data points lie below this line, indicating that the extrapolation is not a good fit. The combined best-fit evolution using all 
the CLASH measurements and our new HUDF12/XDF results is significantly steeper, following (1 + 2) — ll-4±3.l (black solid line). While 
we find an increase by a factor 30 X in the SFRD between z ^ 10 and z ~ 8 from our HUDF12+GOODSS analysis alone, the best-fit trend 
results in a somewhat reduced increase. Nevertheless, this is still a factor lOx in the ~ 170 Myr from z ~ 10 to 2 ~ 8, i.e., a large change 
over a short time period. 



~ 2(7 below the simple extrapolation of the lower redshift 
trends. 

At z 9 our SFRD esti mate is in excellen t agreement 
with the measurement of lEllis et al.l (|2013f) based on a 
photometric redshift selection. At z ~ 10, however, we 
find a significantly lower value, mainly due to our inclu- 
sion of a larger dataset covering a larger area, in which we 
do not find any additional candidate. Nevert heless, our 
measurement is within ~ Icr of the result of lEllis et al.l 



(|2013f) . in particular, after we correct their measurement 
down by a factor two to account for the source that is 
likely the result of a diffraction spike in their z = 9.5 
sample (Figure |9]). 

Also at z ~ 10.7, our upper limit is sign ificantly below 
the z = 10.5 estimate of lEllis et ahl (|2013D . This is likely 
due to the wide redshift range probed by our JHi^q- 
dropout sample, wh ich extends from z ~ 9.5 to z ^ 
11.8 (see Figure jS]). lEllis et al.l H2013D consider a strict 
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boundary of z = 10 to z = 11. 

In order to compute the best-fit evolution of the SFRD 
at z > 8, we combine all z > 8 measurements from 
CLASH with our improved estimates at z ~ 9 — 11, 
together with the p revious z ~ 8 SFRD measurement 
(jOesch et alJl2012bl ). The best-fit evolution falls off very 
rapidly, following oc (1 -I- z)~^^-'^^^-^ . This is shown as 
the black line in Figure [131 As can be seen, this is sig- 
nificantly steeper than the z 4 — 8 trends. By z ~ 10, 
the best-fit evolution is already a factor 5x below the 
lower redshift trend. Therefore, the combined constraint 
on the SFRD evolution from all datasets in the literature 
clearly points to an accelerated evolution at z > 8. 

5. SUMMARY AND CONCLUSIONS 

We have used the new, ultra-deep WFC3/IR data over 
the HUDF field as well as the optical XDF data to pro- 
vide a reliable selection of galaxies at z > 8. The new 
observations from the IIUDF12 program push the depth 
of the Hieo imaging deeper by ~0.2 mag compared to our 
previous data from the HUDF09 survey, and they pro- 
vide additional JHhq imaging. The JHno data are very 
useful for selecting some of the first z 9 and z ~ 11 
galaxy samples using the Lyman Break Technique. Fur- 
thermore , we exten ded our previous search for z 10 
galaxies (jOesch et a l. 2012a) to fainter limits by includ- 
ing this new IIUDF12/XDF data set. Our analysis is the 
most extensive search for z > 8 galaxies to date. 

From our full dataset, we find a total sample of nine 
z > 8 galaxy candidates. Seven of these lie in our z 
9 selection. Contamination is always a central concern 
for high redshift samples and after careful analysis we 
expect that the contamination fraction is small, being 
only about 15 - 20%. 

We found that one of the z ^ 9 sources has a very wide 
photometric redshift likelihood distribution, with an am- 
biguous best-fit at Zphot = 2.2 (using ZEBRA; Zphot = 8.6 
using EAZY) . We therefore exclude this source from the 
subsequent analysis. 

We discover a new z ~ 10 source (at z = 9.8 ± 0.6), 
making it one of the very few galaxies known at this 
very high redshift, just 460 Myr after the Big Bang. 
The highest redshift candidate in our sample is XDFjh- 
39546284 that was previously identified at z ~ 10.4. 
However, these new data (the JHi^q in particular) con- 
strain this galaxy to be at z ~ 11.8, if it is at high red- 
shift. This interpretation is problematic and has led to 
discussion about it being a possible lower-redshift z 2 
object, and so i ts true nature remains quite uncertain 
at this time (see lEllis et al.ll2"0Tl iBouwens et al.ll2012bl: 
iBrammer et al.ll2013f) . 

This sample of z > 8 galaxy candidates proves to be 
very important for setting a number of constraints on 
galaxy build-up at very early times, allowing us to derive 
an estimate of the UV LF at z ^ 9, to improve our 
constraints at z 10, and to set limits at z 11. 

The main result from our analysis is a confirmation 
of our previous finding that the galaxy population, as 
seen down to Mjjv = —17.7 mag, evolves much more 
rapidly at z > 8 than at lower redshift (identified as "ac- 
celerat ed evolution" ; IBouwens et al.|[2011al : lOesch et al.l 
l2012aD . This is seen in (1) the expected number of galax- 
ies when extrapolating the lower redshift trends to z > 8 
(Figure [To)). (2) in the direct constraints on the UV LF 



(Figure [TT]). and (3) in the evolution of the luminosity 
and star-formation rate densities down to our current 
completeness limits (Figures [T2l [13)) . All measurements 
consistently point to accelerated evolution at early times, 
beyond z ^ 8. 

Specifically, if the lower redshift trends of the UV LF 
are extrapolated to z ~ 10, we would have expected to 
see 9 ± 3 candidate sources in our full data set, ^ 5 of 
which only in the HUDF12/XDF data alone. However, 
only one such candidate is found in the HUDF12/XDF 
data, which suggests that the galaxy population evolves 
more rapidly than at lower redshift at 99.5% significance 
(see Section [^?T|) . 

From z ~ 8 to z 9, the luminosity density in star- 
forming galaxies with SFR> O.7M0yr~-^ (i.e. Muv < 
— 17.7) decreases by 0.6 ± 0.2 dex. This decrement is 
fully consi stent with previous estimates from CLASH 
(Bouwcns et al.ll2012a| ) and from the HUDF12 data alone 
( Ellis et aL! ^20T3tl 

The combination of our new measurements of the 
SFRD at z > 8 with all previous estimates from the 
CLASH data (BouwenseLd] l2012at iCoe et all IMsl : 
Zheng et a l. 2012) results in a best-fit evolution which 
is extremely steep, following oc (1 -I- z)~^^-*^^-^ . 

These results on the galaxy evolution at z > 8 con- 
trast with the conclusions drawn by several recent pa- 
pers, who argue that the UV LD evolution at z > 8 is 
cons istent with the lower redshift trends (e.g. lEllis et ahl 
[Ml iCoe et al.ll2013t IZheng et al.ll2012l ). However, the 
small sample sizes of z > 8 galaxies in these individual 
analyses resulted in large uncertainties on the LD and 
SFRD evolution. We show here that once all these mea- 
surements are combined self-consistently, they do indeed 
point to accelerated evolution at z > 8, consistent with 
theoretical expectations. 

Note that the steep fall-off we find in the UV LD at 
z > 8 is not at odds with galaxies driving reionization. 
Our measurements only reach to ~ L* at z ^ 10 (i.e. 
to Muv — —17.7 mag). However, with the steep faint- 
end slopes that are consistently found for z > 4 UV LFs, 
the total luminosity density is completel y dominated by 
galaxies below this threshold (see e.g . [Bouwens et al.l 
[2012(1 [Kuhlen & Faucher-Ciguerdl20ra ). 

With WFC3/IR we are now in a similar situation in 
studying z~9.5 — lOaswe were three years ago with 
NICMOS at z ^ 7. Galaxy samples are still small, and 
the conclusions are uncertain. However, over the next 
few years the z > 8 frontier will be explored more exten- 
sively. In particular, the additional deep field observa- 
tions to be taken as part of the Deep Fields Initiative (a 
large Director's Discretionary program), will significantly 
increase sample sizes and should allow for improved con- 
straints on the z ~ 9 and z ~ 10 LF at Mjjv < —18.5 
mag. This will enable more precise constrains on the 
accelerated evolution that we now see in the galaxy pop- 
ulation from the data over GOODS-South. 
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TABLE 5 

Photometry of 2 > 8 LBG Candidates in the HUDF12/XDF Data 



ID 



RA 



DEC 



^ZEBRA 2B-4zy comments 

pnot pnot 



^160 (yJ)-JHi40 J125 — -f^ieO -^-^140 — -^160 S/N//160 S/Nj/fi4o S/Nji25 



opt 



' 9 YJ-dropouts 



XDFyj-38135540 
XDFyj-39478076 
XDFyj-39216322 
XDFyj-42647049 
XDFyj-40248004 
XDFyj-43456547 
XDFyj-39446317'' 



03:32:38.13 

o .+0.1 



03:32:39.47 



-27:45:54.0 

8 4+°-^ 
°-^-0.2 



27.95 ±0.10 0.8 ±0.1 0.0 ±0.1 -0.2 ±0.1 13.1 16.0 9.8 0.2 

Bouwens UDFy-38125539; McLure HUDF12-3813-5540 {z = 8.3); and in other Y-dropout samples. 



8.1 



+0.3 
-0.6 



-27:48:07.( 
i+0.2 

-0.5 



8.3"' 



28.53 ±0.14 0.8 ±0.2 0.4 ± 0.2 -0.0 ± 0.2 
Bouwens UDFy-39468075; Ellis HUDF12-3947-8076 [z = 



8.7 



8.6) 



03:32:39.21 



8.8 



-0.5 



-27:46:32.2 
5+0.5 

-0.4 



8.9'' 



03:32:42.64 -27:47:04.9 



g 0+05 



03:32:40.24 

o + 0.5 

^'■^-o.s 



n 2+0.5 
^■^-0.6 



29.49 ±0.25 1.1 ±0.5 0.3 ± 0.4 
Ellis HUDF12-3921-6322 (z = 8.8) 

29.15 ±0.21 1.5 ±0.7 0.8 ± 0.6 
Ellis HUDF12-4265-7049 {z = 9.5) 



-27:48:00.4 29.87 ± 0.30 



1.3 ±0.7 



0.3 ±0.6 



0.1 ±0.3 



0.1 ±0.3 



-0.1 ±0.4 



8.9 



+0.6 



03:32:43.45 -27:46:54.7 



8.7- 



+0.6 
-0.5 



8.9 



+0.7 

-0.8 



Faint source. Not in Ellis et al. (2013) sample. 

29.69 ±0.42 1.3 ±0.7 0.5 ± 0.7 -0.1 ±0.4 
Ellis HUDF12-4344-6547 (z = 8.8) 



03:32:39.44 



2.2 



-0.8 



-27:46:31.7 
,+1.0 

-1.5 



29.77 ±0.27 



1.1±0.7 



> 1.0 



-0.3 ±0.5 



5.1 



4.4 



3.5 



3.1 



3.8 



9.7 



4.9 



4.9 



3.3 



3.5 



3.7 



5.0 



3.2 



2.1 



2.2 



1.3 



8.6"* 



Faint source. Very wide p{z), with low best-fit redshift. Not in Ellis et al. (2013) sample. 



-0.6 



-1.4 



-0.5 



0.0 



1.3 



2 ~ 10 J-dropouts 

XDFj-38126243 03:32:38.12 -27:46:24.3 29.87 ± 0.40 > 1.9 1.4 ± 0.9 0.3 ± 0.4 5^8 1^2 -0.6 

9-8^0 6 '^■^io 6 This source was selected as a 2 ~ 10 candidate in the HUDF09 year 1 data, but did not appear in 
final Bouwens et al. (2011) sample due to low S/N in 2nd year data (see Fig|6ll. 



2 '~ 10.7 JH-dropouts 

XDFjh-39546284 03:32:39.54 -27:46:28.4 28.55 ±0.14 - > 2.3 > 2.3 tI 02 To OS 

11.8^04 11-9^0 5 HUDF09 source of Bouwens et al. (2011), Oesch et al. (2012); Ellis HUDF12-3954-6284 (2 = 11.9) 



Note. — S/N are measured in circular apertures of fixed 0''35 diameter 
" Due to the low photometric redshift estimate, we do not include the source XDFyj-39446317 in our analys is of the UV LF at 2 ~ 9. One 
contaminating lower redshift source is expected in our sample due to photometric scatter (see Section 13.5.211 . 
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TABLE 6 

Photometry of Additional Potential z > 8 LBG Candidates not used in this Analysis 



ID 



2, 



RA 

ZEBRA 
phot 



DEC 

EAZY 
^phot 



comments 



{YJ)-JHiio Ji25-J^i60 JHiio-Hieo S/Nhwo S/Nj„^io S/Nji 



^opt 



42126501 



43246481 



43286481 



03:32:42.12 -27:46:50.1 
-+O.I Q 7+0.1 

-0.1 ^-'-0.2 



9.5"* 



28.45 ±0.05 1.7 ±0.2 1.4 ±0.2 0.3 ±0.1 22.2 15.8 

This potential source is completely blended with a foreground galaxy. 



03:32:43.24 -27:46:48.1 28.61 ± 0.17 



0.8 ±0.6 



1.3 ±0.7 



2.5 



+6.1 



2.4 



+ 7.7 



Close to bright, clumpy foreground galaxy. 



03:32:43.28 -27:46:48.1 28.53 ±0.17 



> 1.6 



> 1.9 



0.7 ±0.4 



0.8 ±0.4 



5.4 



6.0 



3.0 



2.5 



10.4; 



-0.5 



10. 



3-0.3 



Close to bright, clumpy foreground galaxy. 



5.2 -0.1 



2.8 -1.9 



0.5 -0.2 



Additional Sources From Ellis et al. (2013) 

UDF12-4106-7304 03:32:41.06 -27:47:30.4 

- - The photometry of this source is significantly affected by a diffraction spike (see Fig. |9ll. 

After subtraction of the expected flux of the spike, the source is only a 2.8(T total NIR detection. 

UDF12-3895-7115 03:32:38.95 -27:47:11.5 30.02 ± 0.30 0.5 ± 0.5 0.1 ± 0.5 0.1 ±0.5 3.5 3.7 3.7 1.2 

0.5^0 5 O.etg 5 Does not satisfy our color selection: (YJ) - JHuo = 0.5 < 0.75. 

Additionally, the best-fit photometric redshift is only 0.5 using our photometry. 



Additional Previous z > 8 Candidates From Bouwens et al. (2012) 



UDFy-37806001 


03:32:37.80 
'•«-0.4 


-27:46:00.1 

70+0.3 


28.39 ±0.12 0.4 ±0.1 
Too blue in {YJ) - JH140 for 


-0.2 ±0.2 
our selection. 


-0.2 ±0.1 


10.1 


13.7 


9.4 


2.7 


UDFy-33446598 


03:32:33.44 
7 7+0.4 
'•'-0.4 


-27:46:59.8 

J 7+0. .5 
' ■ '-0.4 


29.00 ±0.20 0.3 ±0.2 
Too blue in (YJ) - J H 14,0 for 


0.1 ± 0.3 

our selection 


0.1 ± 0.2 


5.6 


6.0 


4.7 


-1.9 



Note. — S/N are measured in circular apertures of fixed 0'.'35 diameter. 



